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INTRODUCTION 


Data analysis In ultrasonic Interrogation often requires quantitative 
methods that provide solution In analytic form. The scattering theory re- 
cently developed within the scheme of an extended version of the equivalent 
inclusion method [1,2] gives such a solution form. Of particular interest 
is that this solution reduces to that of the static solution of the Eshelby 
problem [3], This suggests that many of the research techniques developed 
in the static theories of composite materials and constitutive equations 
can be extended to include the frequency spectra of the ultrasonic waves 
such that the signals received at the probe(s) can be interpreted in terms 
of the microstructural parameters. The far field solutions needed in data 
analysis of ultrasonic evaluation therefore are often complemented by the 
near field solutions obtained in fracture mechanics or theories of composites. 

The determination of the static effective properties of composites 
or poly-crystalline alloys is an area that constitutes extensive body of 
literature. Basically, there are two major theoretical approaches that can 
be used to describe the global material properties which include the causal 
effect introduced by the constituents, i.e. the statistical or deterministic 
approach. Both approaches depend upon the availability of the solution of a 
single Inhomogenelty, be it a crack or an ellipsoidal inclusion material, 
and an averaging theorem. A critique of the approaches and a review can be 
found in [4] and [5], respectively. Recent textbooks on composites are also 
available for references, [6,7]. 

Solution to the scattering of a single embedded inhomogenelty is 
available by different methods that are appropriate at different frequency 
ranges [8]. The methods that offer a solution in an analytic form are: 


the Born approximation [9], the longwave approximation [10-12], the ex- 
tended method of equivalent inclusion [1] , and the extended theory of 
geometric diffraction theory [13,14], 

There are several averaging schemes or theorems that exist in the 
literature for finding the dynamic effective moduli and mass density. Kuster 
and Toksoz [15] considered a sphere of radius "a" consisting of small spheres 
of inclusion material in an elastic matrix. At Rayleigh limit and at low 
concentrations of inclusion material they showed that the scattered displace- 
ments are equivalent to those obtained by a homogeneous sphere with effective 
properties. Berryman [16,17] formulated a self-consistent dynamic theory of 
composites by requiring the scattered displacements to vanish on the average. 
His theory is appropriate for the longwavelength regime. Gaunaurd and 
Uberall [18] studied the resonant scattering from spherical cavities in 
elastic and viscoelastic media. 

Employing the concept of volume averaging for such physical quantities 
as stress, strain and displacements, the moduli and mass density of the ef- 
fective medium are derived by matching the strain energy and the kinetic 
energy of the effective and the physical media. The approach taken is 
similar to those of Christensen and Lo [19], Chou, Nomura and Taya [20] 
and Kohn [21] except the schemes developed in References [19,20] apply 
to the static case, and that in Reference [21] applies to the Rayleigh limit. 
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THEORETICAL CONSIDERATIONS 


Consider the problem of a matrix material, with elastic moduli 
Cjkrs mass density p, that contains a concentration of second-phase 
particles of moduli and mass density p', randomly distributed over 

the material system. A schematic diagram of the system Is as shown in 
Fig. 1. The true composite thus occupies the whole region and possesses 
effective properties and p*. To determine the effective moduli 

and mass density, the following definitions are used: 

<g> = c* <e> 

■^2 = "J C * <£><£> 


•i<p V v> = -^ p*<v><v> 


( 1 ) 

( 2 ) 

(3) 


where a, e and v are the stress, strain and velocity fields, and the 
notation < > denotes the volume average of a field quantity, e.g. 

<F> = (V)"^ / F(x )dV (4) 

V 

in which V represents volume. The left-hand -side and the right -hand- side 
in each of the Equations (1-3) can be shown to be equivalent in a self- 
consistent manner [22,23]. 

Using the elastodynamic solution for a single ellipsoidal inhomo- 
genelty in a linear elastic medium developed in [1], the displacement 
and strain fields inside an inhomogeneity are found, for Incident time- 


harmonic plane wave, to be: 



(5) 


(6) 
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in which the superscripts (m) , (a) and < denote "mis-match, " "applied," 
and "inside the inhomogeneity, " respectively. Employing the volume 
averaging process as described by Eq. (4) and substituting in Eqs. (1-3), 
the effective properties are easily defined as follows: 

p* = p + f Ap D D (7) 

C* = C + f AG A (8) 

where Ap = p'-p, AC = C'-C and f is volume fraction of inclusion material. 
The tensor fields D and A are of ranks two and four, respectively, and 
they are functions of wavenumber, geometric properties, and effective 
and inclusion material properties; 


D . = 6 . - <f .[0] + 4Tr(p'-p*)co^} 

mj mj ' mj ' mj 

A = ■^ M (f) + F (f)\ S., /2p*u)^ + 

mnpq Xmjk,n njk,m / jkpq 


6 6 
mp nq 


(9) 

( 10 ) 


in which the f and F tensors are defined in [1] and is the connecting 

tensor between the elgenstralns and the applied strains, i.e. 


= S 

jk jkpq pq 


( 11 ) 


for the case of uniform eigenstrains and eigenforces. In developing these 
expressions the volume average of the C’-integrals, as defined in References 
[1,2], must be evaluated. Finally, it should be noted that p* and C* are 
physical quantities and hence only the real parts of D D and A should enter 
into consideration in Eqs. (7,8). 
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SPHERICAL INCLUSION MATERIALS 


Let the spherical Inclusion materials of radius "a" be randomly 
distributed over the whole volume of the matrix. If the matrix and the 
inhomogeneities are isotropy, the effective medium is also isotropic. 

It is straight forward to show that 

°mj = ^mj ^^33^°^ (12) 

and 


'jkpq 

^kjpq 

^jkqp 

pqjk 


= s > 

! s «= C 

111 " 

^2222 

3333 

1 


= s 

= s 

= c 

2323 

1313 

1212 

3 

1122 

^1133 

®2233 

= S 


where 

Cl - (C* + C* - 2GC*)/[(C*)^ + C* C* - 2(C*)^] 

= (C*G - C*)/[(C*)2 + C* C* - 2(C*)^] 

C^ = {2 Fj^ 22 1^°^ y*/(p'-p*)} 

^* 1 - G^lll,lt°’ (G+l)D^22,lt0l 

F IT _(X* + 2p*)/G 
G = (A'-X*)/[(X'-X*) + 2(p'-p*)] 

H = X/G 

Following the theory developed in the previous section, the effective 
moduli and mass density are found to be 
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(14) 


p* = p + f Ap 

X* = A + + 2A^j^22^ + (^5) 

p* = p + f 4p(Aj^2i2 '*’ '^1221^ 

K* = K + ft(Aj^^j^^ + 2A^i22^^^ (2/3) [3A^^^22 ^1212 ^1221^^^"^ 

It is clearly seen that the velocities are dispersive. At frequency range 
above that of the Rayleigh limit, this phenomenon is pronounced. From 
Figs. 2-7, the bulk moduli, shear moduli and longitudinal velocities are 
shown as functions of volume concentration of spherical inclusion materials 
for the cases of (1) aluminum spheres in germanium and (2) voids in silicon 
nitride, HS 130. 

As an example of application to detect localized damage by void 
nucleation, let all voids be locally nucleated within a region fi of radius 
R, R>>a, Fig. 8. The effective mdouli of this composite can therefore be 
obtained as before. If void nucleation outside the region fi can be ignored, 
then the scattering of the composite sphere can be easily obtained. Using 
the program developed in [24], the scattering cross section for a composite 
sphere consisted of small voids in titanium is displayed as a function of 
dimensionless wavenumber for different concentration of voids. Fig. 9. 

The scattering cross section, which is essentially proportional to the 
attenuation [25], increases with increasing concentration f. It appears 
that these curves can be used to locate and calibrate porosity in a struc- 
tural component. Effective properties for this material system are pre- 
sented in [26]. 
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CLOSING REMARKS 


The velocity and attenuation of ultrasonic waves in two-phase media 
are studied by using a self-consistent averaging scheme, that require 
the effective medium to possess the same potential and kinetic energy as 
the physical medium. The concept of volume averaging for physical quan- 
tities is employed and the solution depend upon the scattering of a single 
inhomogeneity. The theory is general in nature and can be applied to any 
two-phase material system. Since the scattering of an ellipsoidal inhomo- 
geneity is known, the average theorem presented in this report can be used 
to study the velocity and attenuation of distributed inhomogeneities of 
shapes such as disks, short fibres, etc. The introduction of the orienta- 
tion of these Inhomogeneities besides only their sizes as in the spherical 
geometry will necessarily induce anisotropy in the effective medium. 

Results for randomly distributed spherical inclusions of radius "a" 
are presented. Effective moduli and mass density are found to be dis- 
persive. When the inhomogeneities or voids are nonuniformly distributed, 
attenuation occurs. The case of some localized damage is studied. Since 
it is well known that porosity is directly related to the strength of 
bone [27] and ceramics [28] it appears that the theoretical study of velocity 
and attenuation in two-phase media presents a valuable means for data analysis 
in ultrasonic evaluation of material properties. 
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Fig. 1 


A schematic diagram of an inhomogeneity 
matrix system with effective properties 
and p* 
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Effective bulk modulus vs. concentration: voids 

in silicon nitride 
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Fig. 7 Longitudinal wave velocity vs. concentration: 
voids in silicon nitride. 
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SCATTERING CROSS SECTION 
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Fig. 9 Scattering cross section vs. dimensionless 

wavenumber: distributed voids within a spher- 
ical domain in titanium matrix, a/R=0.01 
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